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We study spherically symmetric black hole solutions with Skyrme hair in the Einstein-Skyrme 
theory with a negative cosmological constant. The dependence of the skyrmion field configuration 
on the cosmological constant is examined. The stability is investigated in detail by solving the 
linearly perturbed equation numerically, ft is shown that there exist linearly stable solutions in the 
branch which represents unstable configuration in the asymptotically flat spacetime. 
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1. INTRODUCTION 

It has been known that the Einstein-Skyrme system possesses black hole solutions as well as regular solutions. 
The black hole solution with B = 1 Skyrme hair was found in Refs. It was shown that there exist two 

• fundamental branches of the solutions and remarkably one of the branches represents stable configuration under linear 
perturbations This is a counter example to the no- hair conjecture for black holes. The no- hair conjecture states 

■ that the only allowed degrees of freedom of stationary black holes are global charges such as mass, angular momentum 
and electric and magnetic charges Q. However, this conjecture concerns only linear field sources and asymptotically 
flat spacetime. In fact, if we consider nonlinear field sources such as non-abelian gauge fields Q and Skyrme fields 
or asymptotically anti-de Sitter spacetime @, black hole solutions with hair can exist. Following the discovery of 
black holes B = 1 Skyrme hair, the B — 2 Skyrme black hole solution with axisymmetry was constructed by us in 
Ref. [T^j- The review of the black hole solutions with Skyrme hair is given in Ref. [TT| ■ All of these solutions are, 

■ however, obtained in the asymptotically flat spacetime. Thus we consider spherically symmetric black hole solutions 
with Skyrme hair in the asymptotically anti-de Sitter (AdS) spacetime with a negative cosmological constant in this 

^ — I ■ paper. 

There has been an increasing interest in the AdS spacetime. Especially the AdS black hole is an interesting 
object from the holographic point of view in the form of AdS/CFT correspondence Maldacena proposed 

the correspondence between the large N limit of conformally invariant SU(N) gauge theory and supergravity (string 
theory) on the product of AdS space with a compact manifold in a low energy limit where the dynamics on the 
brane decouples from the bulk. Witten elaborated on the idea of the correspondence and related the thermodynamics 
of = 4 super Yang-Mills theory in four dimensions to the thermodynamics of Schwarzschild black holes in AdS 
^jTj| space [T3 ■ The study of the AdS black hole in terms of conformal field theory showed that the Hawking-Page phase 
I* ■ transition between AdS space and a AdS black hole corresponds to the confining-deconfining transition in the dual 
. £h ! gauge theory ^4|. It may be thus useful to consider the Skyrme model which is interpreted as QCD in the large N 
^ ■ limit in the asymptotically AdS spacetime to understand the AdS/CFT correspondence. 

Brane world cosmology also indicates that there was a period when spacetime was AdS with a negative cosmological 
constant in the early universe (for example, see |l5Ula .ll7lll8j'). If this scenario is true, AdS black holes should have 
been produced as the primordial black holes. The solutions we obtain in this paper provide a semiclassical framework 
to study the interaction of baryons with the AdS primordial black hole. 

In the context of the Einstein- Yang-Mills (EYM) theory, it was shown that while all black hole solutions are unstable 
for A> [HEl, there exist stable black hole solutions for A < [2l| . We investigate the stability of the AdS Skyrme 
black hole and discuss in detail to see if the cosmological constant change the stability of the AdS black hole skyrmion. 

2. THE EINSTEIN-SKYRME MODEL 

The Einstein-Skyrme system with the cosmological constant A is defined by the action 

S = |d 4 x^j^(i?-2A) + ^ 3 ^tr^ (1) 

where L M = U^d^U and U is an SU(2) chiral field. 

We require that the spacetime recovers the AdS solution at infinity and parameterize the metric as 

ds 2 = ~e 2S(r) C{r)dt 2 + Cir^dr 2 + r 2 {d0 2 + sin 2 6W0 2 ) (2) 
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where 

C( r ) = 1-2^*0-^. (3) 
r 3 

The topology of AdS spacetime is S 1 x R 3 and hence the timelike curves are closed. This can be, however, unwinded 
if we consider the covering spacetime with topology R . For the spacetime to recover the anti-de Sitter spacetime 
asymptotically, solutions satisfy 

m(r) — > const. , 5(r) —* as r — > oo . (4) 
At the horizon r = r^, we have C(r^) = which reads 

mM = ^?~^G ' 5 ^ =Sh (5) 

where Sh is provided by hand as a shooting parameter so as to satisfy the asymptotic boundary condition in Q. 
The B = 1 skyrmion can be obtained by imposing the hedgehog ansatz on the chiral field 

U(f) = cos f(r) + in ■ t sin f(r) . (6) 

Introducing the dimensionless variables 

x = eF n r , n(x) = eF n Gm{r) , A = A/e 2 F 2 (7) 



with 



one gets the energy 



c(z) = i-*a&-^ i (8) 



Es _ ^ [Si [or + Ucr + °m \ ( 9) 
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where we have defined u — x 2 + 8 sin 2 / and v = sin 2 f(x 2 + 2 sin 2 /). 
The covariant topological current is defined by 



\cuf' 2 + -^v J e 5 dx (10) 



BM = -Srr^ tr (u- 1 a u uu-%uu- 1 d t ,u) (11) 

247r^ v~ g 



whose zeroth component corresponds to the baryon number density 



o _ 1 , £an / 
5 --^ e -^5—- ( 12 ) 



r 



We impose the boundary condition on the profile function as 

/(a;) — ► as x — > oo , (13) 

and at the horizon x — Xh, f(%h) = fh is a shooting parameter determined so as to satisfy (|13[) . Then the baryon 
number becomes 



B= f ^B Q d 3 x = -- f sin 2 fdf = J- (2/ h - sin 2/ fc 



(14) 



If the solution is regular and there is no horizon, the B = 1 skyrmion is recovered with //,. = n. Indeed, as shown later 
numerically, the value of fh approaches to n as Xh goes to zero. However, if the solution is a black hole, fh takes the 
value less than 7T, which means that the solution possesses fractional baryonic charge. Topologically, the presence of a 
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spherical event horizon alters the spacetime topology to be R 2 x S 2 and unwinds the skyrmion. Hence the skyrmion 
absorbed partially by the black hole is topologically trivial. 

The Einstein equations with a cosmological constant A takes the form 



(15) 



which reads 



Goo = 87rGT o - A.g o 
G u = 8nGT n - Ag u 



l-C-C'x 



Cuf 
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2v 



-1 + C 



(e 2S C)' 



-..26 



Ax 2 



2v 



x = — [ Cuf 2 -—)-Ax 



(16) 
(17) 



where we have defined the coupling constant a — AirGF 2 . Consequently following two equations are obtained for the 
gravitational fields 



- {Cx y + i = ^(cuf 2 + ^ 



Ax 2 



The variation of the static energy l|l(J|) with respect to the profile f(x) leads to the field equation for matter 



e°Cu 



-{e 6 Cu)'f' + 4GT + 1 + 



4 sin 2 / 



e 6 sin 2/ 



Thus the coupled field equations to be solved are given by 



V = 
f" = 



a I „ ,, 2 2v 



Cuf 2 + 



? s Cu 



-{e d Cu)'f + [ACf 2 + l + 



4 sin / 



e s sin 2/ 



(18) 

(19) 

(20) 
(21) 

(22) 



In order to determine the boundary conditions on the regular event horizon, let us expand the fields around the 
horizon Xh 



2 



Ax 



G 



* + m(x - x h ) + 0((x - x h ) 2 ) 



f = f h + f(x - x h ) + 0((x - x h ) 2 ) 
S = S h + Si(x - x h ) + 0((x ~ x h ) 2 ) . 



Inserting them into the field equations l]2U p -l]22 p . one obtains 



h 



1 + 



2 sin 2 /?, 



sin f h 



4 sin 2 f h 



x h {x 2 h + 8 sin 2 f h )(l - 2/n - Aa: 2 ) 



sin 2 f h 



Si = — (x 2 + 8 sin 2 f h ) ft. 

4x/i 



(23) 

(24) 
(25) 



(26) 

(27) 
(28) 



The profile function numerically computed for several values of |A| are shown in Fig.^ The horizon radius and the 
coupling constant are fixed with Xh = 0.1 and a = 0.02. There are two branches of solutions for each value of the 
cosmological constant as well as the coupling constant as was seen in the asymptotically flat spacetime case. We define 
the solution with larger (smaller) values of fh as upper (lower)-branch. It is shown that the skyrmion shrinks as |A| 
becomes large in the upper branch. In Refs. j^Q, one can see that the skyrmion shrinks as the coupling constant 
becomes large in the upper branch. Therefore increasing the value of |A| gives similar effects on the skyrmion as 
increasing the value of the gravitational constant. On the other hand, in the lower branch, the change in size is much 
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smaller and is almost unrecognizable. In the asymptotically flat spacetime, the solution in the lower branch expands 
in size as the coupling constant increases. Fig. [2] shows the value of fh as a function of x%. It is observed that fh 
continuously approaches to n as Xh goes to zero, recovering the B — 1 regular skyrmion solution at Xh = 0. 

The black hole mass-horizon radius relation is shown in Fig. [3] It has been shown that the area of an AdS black 
hole with scalar hair proportional to the entropy |22l |23| . Thus we infer that this relation still holds for AdS black 
holes with Skyrme hair as 

s =5l = ^)' <29 » 

One can see that the upper and lower branch correspond to the high- and low-entropy branch respectively. The 
cosmological constant reduces the entropy of the black hole. The reduction of the entropy is also seen when the 
coupling constant increases. 

Fig. 0] shows the parameter fh as a function of |A| for a = 0.0, 0.02, 0.04 with Xh = 0.1 fixed. The value of fh is 
directly related to the baryon number as can be seen from Eq. I|14[l . Thus, in the upper branch, the baryon number 
becomes smaller as |A| becomes larger, which represents the baryon more absorbed by the black hole. In the lower 
branch, the baryon number slightly increases as |A| becomes large. This result also shows that the cosmological 
constant gives a similar effect on the skyrmion as the coupling constant. 

We found the maximum value of |A| above which there exists no solution for each value of the coupling constant. 
In Fig. |SJ the maximum value of |A| is shown as a function of a. The maximum value decreases as a increases, 
and at a = 0.126 it becomes zero. Thus at a = 0.126, the asymptotically AdS solution does not exist and only the 
asymptotically flat solution exists. 



3. LINEAR STABILITY ANALYSIS 



In this section we shall examine the linear stability of the black hole solutions described in the previous section. 
Let us consider the time-dependent small fluctuation around the static classical solutions fo, So and fio 



From the time-dependent Einstein-Skyrme action 



f(r,t) = fo(r) + h(r,t) 
5(r,t) = S (r) + S^t) 



(- 



5 C 



f 2 + Cf' 2 )u + v 



e s dx, 



one obtains the time-dependent field equation as 



1 ( 1 



(e°Cuf>y + - \ -j(jf ~ e°Cf 

where we have defined Uf = Su/5f and Vf = Sv/Sf. 
The time-dependent Einstein equations are derived as 



e s v f 1 •• 
u f - —T- = — u -f 



e s C 



Goo — 87rGToo 
Gn =8vrGTii -1 + C 
which reads the following two equations 

S' = 



l-C - C'x = — 
4 



1 



(e 2S C )' 

„2S 



e 2& C 
a 



2v 



U/ 2 + G/' 2 



,25 C 



2v 



1 



Ax \ e 25 C 2 



rf 2 + f 2 



-{Cx)' + 1 



2x 2 



v + CS'x . 



(30) 

(31) 
(32) 



(33) 



(34) 



(35) 
(36) 

(37) 
(38) 
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Substituting Eqs. P0 fl -ip2 )l into Eqs. (EH and gives the linearized equations 



2x 2 



Eq. I|39|l and the classical field equation derived from Eq. I|34[) which is 

are inserted into Eq. H4()(l and resultantly one gets 

- (e^Cxx)' = |(e 5 °aruo/o/i)' 
which can be integrated immediately to obtain 

Ci = -—C u f^ 1 . 

Similarly let us linearize the field equation l|34|> . Using Eqs. (|39[1 . (|41|l and l|43|) . one arrives at 



(e 5o C u /0'-C/o/i 



(39) 
(40) 

(41) 
(42) 
(43) 

(44) 



where 



1 Jon o, f' 2 a_ e "^//o 
'2 e C u //o/0 + ^2 



Setting /i = ^(i)e mt /V«o , we derive from Eq. (glj as 



(e 5 °C O' + 



w 



e a °C 



Let us introduce the tortoise coordinate x* such that 

dx* 
dx 



1 



with — co < x* < +oo. Eq. (|46|) is then reduced to the Strum-Liouvillc equation 

d 2 t 



dx* 2 



where 



Uo = e So C 



2Juq~ 



U £ = u'Z 



e So C 



-Uo 



(45) 



(46) 



(47) 



(48) 



(49) 



The solution is linearly stable if there exist no negative eigenvalues since imaginary u> represent exponentially growing 
modes. Unfortunately the potential Uq has a complicated form and we are unable to discuss the stability analytically. 
Thus we solve the wave equation (I46f) numerically under the boundary conditions £ — ► as x ^ Xfi and x — ► oo which 
ensure the norm of the wave function to be finite. We show the wave function in Fig.|H|for |A| = 1 with a — 0.02 and 
Xh = 1 fixed. The ground state corresponds to the wave function with no node. The first excited state corresponds 
to the wave function with one node. Remarkably both the lower and upper branch solution have ground states with 
oj 2 positive and hence they are stable. Fig. [3 shows the wave function of the lower branch for |A| = 0.5. The ground 
state has a negative eigenvalue (lo 2 = —0.273) and thus the solution is unstable. We have found the critical value of 
|A| at which the stability of the lower branch solution changes and shown in Fig. [S] as a dotted line. For |A| = 0, the 
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lower branch solution has one negative mode and for u 2 > 0, the mode becomes continuous (see Fig. [SJ. Hence it is 
unstable for all values of the parameters a and Xh , which confirms the previous results in Ref. [^, . 

The results we have obtained indicate that the cosmological constant stabilizes the black hole skyrmion. This 
is consistent with the result of the EYM black hole where it was shown that the EYM black hole is unstable for 
A > 0113, but the stable EYM black hole exists for A < HJ. 

Another important feature for the black hole skyrmion with A < is that only discrete modes exist in both of 
the branches. It is because the potential behaves asymptotically as x 2 , meaning all the eigenvalues are discretized 
analogous to the harmonic oscillator eigenvalues. 



4. CONCLUSIONS 



We have studied the black hole skyrmion with a negative cosmological constant. There exists two fundamental 
branches of the solutions; the upper- and lower-branch corresponding to the high- and low-entropy branch respectively. 
The increase in the absolute value of the cosmological constant |A| gives similar effects on the skyrmion as increasing 
effectively the value of the gravitational constant fixing the pion decay constant to be the experimental value. The 
skyrmion shrinks and the baryon number is more absorbed by the black hole as |A| increases. There is the maximum 
value of |A| above which no solution exists for each value of the coupling constant. In particular, we observed that at 
a = 0.126, no asymptotically AdS solution exists and only the asymptotically flat solution exists. 

The linear stability was examined in detail by solving the linear perturbed wave equation numerically. In the 
asymptotically flat case, the upper branch is stable and the lower branch is unstable. However, in the AdS case, 
surprisingly there exist stable solutions even in the lower branch depending on the value of |A|. The observation that 
the negative cosmological constant stabilizes the black hole skyrmion solution is consistent with the result for the 
EYM black hole where there exist stable solutions on ly w hen A < 0, Efl l2ll | . This implies that the catastrophe 
theory for the stability analysis of hairy black holes [2J| may not be applicable to the black hole solution with a 
negative cosmological constant. 

The solutions we obtained in this paper provide a semiclassical framework to study the interaction of baryons with 
the AdS prim ordial black hole. If the AdS primordial black holes were created in the early universe as indicated in 
Refs. |l5l lift ll7L llSj. it would have induced baryon decay. The Einstein-Skyrme theory should be useful as a simple 
framework to study such decay process. 

Since the Skyrme model corresponds to QCD in the large N limit, it may be worth understanding the Skyrme 
model in the context of the string theory. With the assumption that the brane is AdS, the extension of recently 
discovered brane-skyrmions to AdS brane-skyrmions would be one of the possibilities |25| . 
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FIG. 1: The profile function / as a function of the radial coordinate x for |A| = 0.0, 0.5, 1.0 with xh = 0.1 and a = 0.02 fixed. 
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FIG. 2: The value of the profile at the horizon fh as a function of xu with a — 0.02 and |A| = 1.0 fixed. 
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FIG. 5: The solid line shows the a dependence of the maximum value of |A| above which there exists no black hole solution. 
The dotted line shows the a dependence of the value of |A| above which the lower branch solution changes its stability to 
become stable. 
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FIG. 6: Wave function £ for the upper and lower branch with |A| = 1 and a = 0.02. 
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FIG. 7: Wave function £ for the lower branch with |A| = 0.5 and a = 0.02. 
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FIG. 8: Wave function £ for the lower branch with |A| = 0.0 and a — 0.02. 



